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ABSTRACT 

Context. The mechanisms governing the star formation rate in spiral galaxies are not yet clear. The nearby, almost face-on, 
and interacting galaxy M51 offers an excellent opportunity to study at high spatial resolutions the local star formation laws. 
Aims. In this first paper, we investigate the correlation of H2, Hi, and total gas surface densities with the star forming activity, 
derived from the radio continuum (RC), along radial averages out to radii of 12kpc. 

Methods. We have created a complete map of M51 in 12 CO 2-1 at a resolution of 450 kpc using HERA at the IRAM-30m 
telescope. These data are combined with maps of Hi and the radio-continuum at 20cm wavelength. The latter is used 
to estimate the star formation rate (SFR), thus allowing to study the star formation efficiency and the local Schmidt law 
Esfr oc E" as . The velocity dispersion from CO is used to study the critical surface density and the gravitational stability of 
the disk. 

Results. The total mass of molecular material derived from the integrated 12 CO 2-1 intensities is 2 10 9 M©. The 3cr detection 
limit corresponds to a mass of 1.7 10 s Mg. The global star formation rate is 2.56 M©yr -1 and the global gas depletion time 
is 0.8 Gyr. Hi and RC emission are found to peak on the concave, downstream side of the outer south-western CO arm, 
outside the corotation radius. The total gas surface density E gas drops by a factor of ~ 20 from 7OM0pc -2 at the center 
to 3Mqpc -2 in the outskirts at radii of 12 kpc. The fraction of atomic gas gradually increases with radius. The ratio of 
Hi over H2 surface densities, Ehi/Eh2, increases from ~ 0.1 near the center to ~ 20 in the outskirts without following 
■ a simple power-law. Ehi starts to exceed Eh2 at a radius of ~ 4 kpc. The star formation rate per unit area drops from 
~ 400 Mq pc -2 Gyr -1 in the starburst center to ~ 2 M© pc -2 Gyr -1 in the outskirts. The gas depletion time varies between 
0.1 Gyr in the center and 1 Gyr in the outskirts, and is shorter than in other non-interacting normal galaxies. Neither the Hi 
surface densities nor the H2 surface densities show a simple power-law dependence on the star formation rate per unit area. 
In contrast, E gas and Esfr are well characterized by a local Schmidt law with a power-law index of n = 1.4 ± 0.6. The index 
equals the global Schmidt law derived from disk-averaged values of E gas and Esfr of large samples of normal and starburst 
galaxies. The critical gas velocity dispersions needed to stabilize the gas against gravitational collapse in the differentially 
rotating disk of M51 using the Toomre criterion, vary with radius between 1.7 and 6.8 kms -1 . Observed radially averaged 
dispersions derived from the CO data vary between 28 kms -1 in the center and ~8 kms -1 at radii of 7 to 9 kpc. They exceed the 
critical dispersions by factors Qgas of 1 to 5. We speculate that the gravitational potential of stars leads to a critically stable disk. 
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1. Introduction 

M51 is an interacting, grand-design spiral galaxy at a dis- 
tance of only 8.4 Mpc seen nearly face-on (Tabled]). It is 
rich in molecular gas, most of which is found in the two 
very prominent spiral arm s which are presumably caused 
by the tidal interaction (|Tullvl Il974ri iHoward fc Bvrdl 
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1990). Kinematic studies reveal unusually large streaming 
motions implying a strong density wa ve and the presence 
of galactic shocks jAalto et aljfl999^ . The [C 11] map of 
iNikola et alJ l)200lj) shows two secondary lobes of emission 
to the northeast and southwest near the corotation radius 
of the density wave pattern, presumably due to cloud- 
cloud collisions, stimulating star formation. ICalzetti et all 
(|2005^l combine Spitzer infrared data of the Infrared Array 
Camera (IRAC) and the Multi Imaging Photometer for 
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Fig. 1. Map of 12 C0 2-1 integrated intensities [Kkms -1 ] 
showing M51, i.e. NGC5194 and its companion galaxy 
NGC5195 in the north-east. The image has a resolution 
of 11" a nd is constructed from a masked moment cal- 
culation ijAdler et al.lll992|) to minimize the noise contri- 
bution from emision-free channels when integrating over 
the full kinematic extent of M51, 350 kms -1 < vi sr < 
600 kms -1 . Color coded intensities are on a logarith- 
mic scale. The cross marks the 0/0 center position at 
a = 13^29™52?7, S = 47°11'43" (eq.2000). The noise level 
is la = 0.65 Kkms" 1 . Contours range in units of la from 
1, 3, 6, 16 to 56 in steps of 10. The peak intensity is 
47.3 Kkms -1 . All intensities are on the T£ scale. 



Spitzer (MIPS) with UV data from the Galaxy Evolution 
Explorer (GALEX). These data are used to discuss the 
various tracers of star formation in M51. Enhanced star 
formation activity at the two lo bes seen in [C ill is proba- 
bly triggered by the interaction jNikola et al.ll200lh. This 
is also indicated by disk simulations of lToomre fc Toomrel 
lll972l) which show tidal tails emerging from the disk at the 
positions of enhanced star formation. VLA observations of 
the H I line show a prominent tail outside the main disk ex- 
tending to the south and east, and c overing more than 25', 
i.e. 100 kpc, in projec te d distance llRots et al.lll990|) (cf. 
iTilanus fc Alien! Il99l|) . iMeiierink et al.1 l|2005[) find that 
a major fraction of cold dust emission at 850 pm stems 
from an extended exponential disk with a scale height of 
5.45 kpc, possibly also tracing total gas column densities. 
The kinematics of Ha e mission of M51 were recently stud- 
ied by the S INGS team foaigle et alJl2006^ . improving on 
the work of lTilanus fc Alle nl lll99ll) . 

Its proximity and inclination make M51 an ideal tar- 
get to study the efficiency at which it forms stars from the 
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Fig. 2. VLA map of integrated Hi intensities [Jy/beam] 
at 13" resolution l|Rots et al.lll990|) in colors. Contours 
show integrated 12 CO 2-1 intensities (cf. Fig.Q. 



gas under the influence of the interaction, and to compare 
with the star formation laws commonly found in spiral 
galaxies. Spiral galaxies often follow two empirical laws: 
as star formation is fueled by the interstellar gas reservoir, 
the star formation rate (SFR) is proportional to a power of 
the total gas surface density: Esfr oc £™ as . This appears 
to hold globally when E gas and Esfr are averaged over 
the entire star-f orming region of a galaxy l|Schrnirldll 95fll 
lKennicuttlll998l and others) with n = 1.4. This may be 
understood assuming that SgpR is proportional to the gas 
surface density over the free-fall time rg- and that the free- 
fall time and the gas surface density are simple functions 
of the average gas density p: oc p^ 1 ^ 2 and £ gas oc p. 
This simple model l eads directly to t he Schmidt law with 
an index of 1.5 fe.g. lElmegreenlfl994|) . More recent studies 
show that the global Schmidt law also appears to hold lo- 



cally in individual spiral galaxies w hen studying radially 
averaged profiles of E gas and Xsfr 

llWong fc Rli t7i l2002t 



iBoissier et al.ll2003|l . The second law is that the total gas 
surface density i s of the order o f a critical gas surface den- 
sity give n bv thelToom rd lll96 4|) criterion for gravit ational 
stability l|Martin fc KennicutdboOlt iLi et al.ll2006|) . 

The relative importance of the different components 
of the interstellar medium was addressed bv lNikola et alJ 
d200ll) when trying to explain the mapped [Cn] emis- 
sion. In a recent study comprising all major cooling 
lines of the molecular gas, i.e. those of [Cn], [Oi], CO, 
and [C i] , at selected spiral a rm positions and the center 
of M51. iKramer et all l)2005|) conclude that the bulk of 
the emission stems from clumpy photon dominated re- 
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Fig. 3. Map of radio continuum intensities at 20 cm 

(VLA C+D, 1 5" resolution) in units of Jy/beam 

(jPatrikeev et all l200fil) . Contours show integrated 12 CO 
2-1 intensities (cf. Fig.QJ. 

gions. Only about 15-30% of the [Cn] emission stems 
from dense Hn regions. The [Ci] and CO emission of 
the ce nt er region have also b e en studied bv llsrael et all 



ll200flh: llsrael Raasl J2002I) : ICerin fc Phillinsl l|2000l) . 
ISchinnerer et alJ l|2004h mapped two distinct regions in 
the spiral arms of M51 in HCN, HCO + and other tracers 
of the chemistry. 

iRvdbeck et alJ l)2004() presented a large CO 1-0 map 
of M51 obtained at the Onsala-20m telescope. The in- 
ner region of M51 were previously mapped in CO 1— 
with the F CRAO-14m, IRAM-30m and NRO-45m 
telescopes bv IScoville fc Yound lll983tk lLord fc Yound 
Jl990tl:lGarcia-Burillo et alJ dl993rJk iNakai et all l|l994|) ; 
iKuno et all <|l995h : iKuno fc Nakail l|l997l) . These single- 
dish observations reveal the large-scale emission for which 
intcrfcrometric observations are not sensitive. OVRO was 
use d to create CO 1-0 maps bylRand fc Kulkarnil l)l990j) 
and lAalto et~aXl (|l999() . iHelfer et etl.1 (|2003^ used BIM A in 
combination with NRAO-12m single-dish data to map the 
inner part of M51 in CO 1-0. 

Here, we present a complete IRAM-30m CO 2- 
of M51, improving on the previous data of 



1 map 01 Mol, i mproving: 
iGarcia-Burillo et alJ l|l993bla|) 



2. Observations 

Observations of the 12 CO 2-1 emission from M51 were 
conducted with the IRAM-30m telescope in February 2005 



Table 1. Basic properties of M51. References: a RC3 cat- 
alogue |de^^^^^^^^n_|^p3)> b iFeldmeier et alJ dl997f) . 
note that iTakats fc Vinkol (l2006t) recently reported a distance 
of only 7.1±1.2Mpc, c |Tuii3 |l974a|) ■ 



RA(2000) 

DEC(2000) 

Type 

Distance [Mpc] 
11" correspond to 
Heliocentric velocity [kms~ 
Position Angle [deg] 
Inclination [deg] 



M51 

13:29:52.7 
47:11:43 
SA(s)bc pec a 
8.4 b 
448 pc 
463 a 
170 c 
20 c 



using th e 18 element focal pla ne heterodyne receiver array 
HERA l|Schuster et al.ll2004[) together with the WILMA 
autocorrelator backend. WILMA has a channel spacing 
of 2.6 kms" 1 (2 MHz) and a bandwidth of 1200 kms" 1 
(930 MHz). Pixel 2 of the second HERA polarization 
showed excess noise and was ignored from further anal- 
ysis. 

Observations were conducted in position switched on- 
the-fly (OTF) mode scanning M51 in right ascension. 
Sampling was 6" in RA. HERA was rotated by 18.5° to 
obtain a s pacing of 7.6" in dec lination between adjacent 
scan lines ^Schuster et alJEooih . This corresponds to near 
Nyquist sampling for a half power beamwidth (HPBW) 
of 11". The resulting map has a size of 11' x 11'. Figure^ 
shows the 7' x 10' sub-region where CO was detected. 
An emission-free reference position was selected at offsets 
(10', 0'). The mean baseline rms is 18 mK at 5kms _1 ve- 
locity resolution on the Tj£ scale. To correct to first order 
for the telescope efficiencies, i.e. to go from the T£ scale to 
the T m b scale, we simply multiplied the antenna temper- 
ature data with the ratio of forward efficiency F e g = 0.91 
over main beam efficiency -B e ff = 0.52. These numbers 
show that the mapped spatial structure of M51 is to some 
extent sm eared out by the e rror beams of the IRAM-30m 
telescope ijGreve et al.lll998|) . All data reduction was done 
using the GILDAS 1 software package supported at IRAM. 



3. Data 

3.1. Molecular gas distribution 

The HERA map of 12 CO 2-1 (Fig.QJ is the first CO map 
of M51 encompassing the companion galaxy as well as 
the south-western arm out to a radius of ~ 12 kpc in a 
homogeneously sampled data set at linear scales of down 
to 450 pc. 

The emission detected with the 30m telescope traces 
the well known two-armed spiral pattern out to the com- 
panion galaxy NGC 5195, which shows up brightly in the 
north-east at ~ 10.5 kpc radial distance, and out to the 



1 |http : //www . iram . f r/ IRAMFR/GILDAS 
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south-western tip of the second arm at the opposite side 
of M51. The outer parts of the two arms in the west and 
in the east appear more fragmented than the inner parts. 
The western arm especially is almost unresolved. Inter- 
arm emission is detected above the 3cr level out to radii 
of about 6kpc. Several spoke- like structures connect the 
spiral arms radially. 

To estimate the total H2 column densities from the in- 
tegrated CO 2-1 intensities , we assume a 2-1/1-0 inten - 
sity ratio of 0.8 as found bvlGarcia-Burillo et alJ 111993a) . 
iGarcia-Burillo et all l|l993b|) and iGuelin et all ljl995|) in- 
dependently derived the CO-to-H 2 conversion factor X, 
for M51 from CO 1-0 and dust continuum data. They 
find that it is a factor 4 smal l er than the Milky Way 
value X M w (jStrong et alJll988t Istrong fe Mattoxlll996|) 
and const ant with radius . Simil ar X-factors for M51 were 
found bv iNakai fc Kunol <) 19951) . using CO data and vi- 
sual e xtinctions towards H11 regions, and bv lAdler et alJ 
l|l992() using BIMA CO 1-0 data assuming that the GMCs 
are in virial equilibrium. 

For the total H2 column density per beam in M51, we 
use: N(U 2 ) = 0.25 X M w (1/0.8) / T mb (CO(2 - l))dv with 
I MW = 2.3 10 20 cm _2 (Kkms- 1 )- 1 . 

The X-factor may be a function of the metallicity in 
spir al galaxies as has be en suggested by several authors 
(e.g. lArimoto et al~lll99(f) . However, other factors like the 
radiation field or the cosmic-ray rate also have a strong im - 
pact on the CO-to-H 2 conversion factor (|Bell et alJl200q) . 
The metallicities of M51 have been found to be slightly su- 
persolar showing only a shallow drop with radius by only 
-0.02dexkpc~ 1 (|Bresolin et alJl2004j) : 12 + log(0/H) = 
8.72(±0.09) - 0.28(±0.14)i?/i? with R = 5.4'. The al- 
most constant metallicity appears to be consistent with a 
constant X-factor. 

The 3<T limit with resolutions of 11" and 5kms _1 cor- 
responds to a mass of 1.710 5 M Q . The spatial resolution 
of 450 pc does not allow to detect individual GMCs if 
their typical size is ~ 50 pc. The individual "clumps" de- 
lineating the spiral arms like beads on a string (Fig.GJ 
have been labeled giant molecular associations (GMAs). 
Th ese may be bound clu sters of GMCs as suggested 
bv lRand fc Kulkarnil l|l990|) or random superpositions o f 
molecular clouds and GMCs (|Garcia-Burillo et alll993a| . 

The total mass of the molecular gas of M51 derived 
from the CO 2-1 data set is 1.9410 9 M Q . This value 
agrees well with th e total mass of 1.610 9 M Q derived by 
iHelfer et all (|2003j) from NRAO 12m CO 1-0 data when 
using the same CO-to-H2 conversion factor and distance. 

3.2. H 1 and radio continuum 

The large-scale distribution of the 21 cm line of ato mic 
hydrogen in M51 was analyzed bv iRots et all l)l990l) us- 
ing the VLA. Note that the total flux is in reason - 
able agreement with single-dish observations l)Rotslll980l) . 
The Hi emission at 13" resolution (Fig.[5J) is weak in 
the inner region while the outer CO arms are clearly 



delineated in Hi. There is patchy Hi emission in the 
south, near (0,-250"), where only very little CO is found. 
In the inner a rms of M51, inside the cor otation radius 
of rsj 7.4 kpc llGarcia- Burillo et aflll993al). Hi and Ha 
emission llRand fc Kulkarnil Il990t iTilanus fc Allen] Il99li 
IScoville et al.ll200ll) as we ll as young star cluster com- 
plexes IjBastian et al.ll2o"o5j) are seen slightly towards the 
convex side, i.e. downstream relative to the CO emission, 
suggesting that they arise when GMCs are destroyed by 
short-lived OB stars. Figure|21 shows that the Hi clouds 
tracing the south-western spiral arm, peak on the concave 
side of the CO arm. This is outside the corotation radius 
and thus again downstream, consistent with the above in- 
terpretation. 

Figure^ shows a map of the radio continuum (RC) 
at 20 cm (1.4 GHz) taken with the VLA at 15" resolu- 
tion (|Patrikeev et al.ll2006|) . Further below, we will use 
this map to derive the star formation rate per unit area. 
The radio continuum map shows strong emission in the in- 
ner arms, which slowly drops off towards the outer arms. 
There is patchy 20 cm emission in the south where there is 
H 1 but only very little CO. Again similar to H 1, the 20 cm 
emission tracing the south-western spiral arm peaks on 
the concave, downstream side of the CO arm. The spiral 
arm near (—40", 125") shows an interesting discontinuity 
and local lack of radio continuum, Hi, and CO emission. 
IPatrikeev et alJ l|200q) discuss systematic offsets between 
the spiral arms traced by 6 cm radio continuum emission, 
ISOCAM data at 15 /im, and BIMA CO 1-0 emission over 
scales of several kpc. 
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Fig. 4. Radial distributions of surface densities of H2, Hi, 
the total gas, and the star formation rate per unit area 
in M51. The lower box shows the ratios E gas /EsFR and 
Eh 2 /Esfr- 
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Fig. 5. Radial distribution of the ratio of atomic J2m ^° 
molecular surface densitie s for M51 (thick line), the 
six spirals of the sa mple of lWong fc Blitd <|200^ . and M33 
<)Hever et al.ll2004T) . For comparison, we show two power- 
laws, one where the ratio scales with R 15 and one where 
it scales with R. 



J^IO 




10 - 



1 10, 
S [M pc^] 

Fig. 6. Radially averaged star formation rate per unit 
area, X)sfr> versus surface densities of the total gas X) g as' 
and of H2 and Hi only. The solid green line is the lo- 
cal Schmidt-law found in M51. The d ashed-dotted blac k 
line is the global Schmidt-law found bv lKennicuttl l|l998h . 
Drawn blue lines represent lines of constant gas depletion 
time or star formation efficiency. 



4. Results 

We investigate the correlation of H2 and H 1 surface den- 
sities with star forming activity along radial averages, in- 
cluding the outskirts of the disk where Hi dominates. 

Various tracers of star formation such as Ha, UV, the 
far infrared continuum, and the radio continuum are fre- 
quently used. Ha emission is directl y linked to OB associ- 
ations but subject to extinction. See lBastian et all l)2005|) 
for a recent study of the local star formation rate in M51, 
derived from Ha images. Optically thin far-infrared (FIR) 
emission is usually taken as the m ost direct indicator o f 
star formation. In a recent study, ICalzetti et alJ l)2005[) 
used MlPS/Spitzer 24 /xm maps to study the star forma- 
tion rate. Some questions remain however in regions of 
diffuse FIR emission from the atomic ISM. In such regions 
the dominant contribution to the FIR flux may be due to 
radiation h eating of dust by the normal interstellar radi- 
ation field lICox fc Mezgerlll989h and thus not trace star 
formation. 

The radio continuum has widely been accepted as 
an alternative measure for star formation activity, a fact 
which is underlined by the very strong FIR/RC correla- 
tion i n many different galaxies (see the review bv lCondonl 
Il992fl . The FIR/RC correlation holds also on local scales 
of a few hundred parsecs comparable to the resolution of 



our data set. See lMurohv et al who recently stud- 

ied this relat ion in M51. Here, we use a new large scale 
20 cm map of lPatrikeev et al"l 1 2006t) at 15" resolution and 
assume that the 20 cm radiation can be used as a direct 
indicator for the star formation rate. We use the FIR/RC 
correlation and a subsequent SFR/FIR conversion to de- 
rive the star formation rate. We then compare the star 
formation rate with the H2 and H 1 surface densities in 
order to check for star formation thresholds. We further 
use the measured CO 2-1 line widths to study the gravi- 
tational stability of the gas disk. 

4.1. Radial distribution of the gas 

Total H2 column densities are derived from the integrated 
CO 2-1 intensities as described in Sec. 13. II We derive 
the face-on surface density via = 2mnN(H2) cosi 

with the atomic hydrogen mass mu and the inclina- 
tion angle i of M51 (Tabled). The total Hi col- 
umn density is derived assuming optically thin emission: 
iV(HI) = 1.82 10 18 /T mb du/(Kkms _1 )cm- 2 . The corre- 
sponding surface density is: X)hi = m H -/V(HI) cosi. The 
total molecular mass of M51 is 1.94 10 9 M Q (SecEHJ and 
the global ratio of Hi over H2 mass is 1.36. The molec- 
ular gas content of M51 is similar to the total molecula r 
mass of the Milky Way, 1.3 10 9 M Q l|Misiriotis et all2006|) . 
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However, the Milky Way has a larger fraction of H I mass. 
The ratio of Hi over H 2 mass in the Milky Way is 6.3. 

Figure^] shows the radial profiles of the surface den- 
sities which were created by averaging in elliptical annuli 
spaced by 10", including points of no detection, and cen- 
tered at the 0/0 position. The annuli were assumed to be 
circular rings viewed at an inclination angle i and with a 
line of nodes rotated from north to east by the position 
angle PA as listed in Tabled 

Total gas surface densities are given by yj„ as = 
1.36 (y^u + y^jji) which includes helium. The radially av- 
eraged yj gas drops by a factor of ~ 20 from 70 M Q pc~ 2 in 
the center to3M pc -2 at radii of 12.2 kpc (300") (Fig.©. 

Figure|5] shows the variation of Ehi / Eh 2 w ^h radius 
in M51. Standard deviations of the mean yj HI and 
are calculated from the rms values along t he annular av- 
erages : a = rms (Apix/n.boam) -0 ' 5 ; following lWong fc Blit d 



( 2002) with the number of pixels per annulus iV p i x and the 
number of pixels per beam riboam- The fraction of atomic 
gas gradually increases with radius. In the inner regions 
of M51, up to ~ 6 kpc, the ratio of HI to H2 surface den- 
sity scales roughly with R 15 . As the radial HI profile in 
M51 is roughly constant, varying only by a factor of ~ 5, 
the power-law dependence of the ratio of molecular and 
atomic gas is largely due to the decrease of the molecular 
gas surface density which drops by more than 2 orders of 
magnitude. The Hi gas surface density starts to exceed 
the molecular surface density at radii greater than 4 kpc. 
The outer regions show deviate from any simple power- 
law. The peak E HI I Eh 2 fraction reaches values of 20 in 
the outer areas of our map. 

We also plot the fraction of atomi c gas for the seven 



spirals studied bvlWong & Blitd" l|2002l) and for M33 stud- 
ied by iHever et alJ l)2004l) . In the center regions, the 
Ehi / Eh? ra ti° 01 the nine galaxies varies strongly be- 
tween 0.03 2 in NGC4321 and slightly more than 1 in M33. 
With the exception of M33, the inner regions are clearly 
dominated by molecular material. For the galaxies of the 
Wong sample, the ratio scales roughly with R 1 - 5 for radii 
upto ~ 6 kpc. In M33, the atomic gas dominates the total 
gas surface density for the entire disk, in strong contrast 
to the other 8 spirals (Fig-EJ. The slope is much more 
shallow: the fractio n of atomic gas scales with R 6 only 
l|Hever et alJbOO^ - At — 7 kpc, the fraction equals the 
fraction found at that distance in M51. 

Next, we compare the total gas surface densities with 
the star formation rate derived from the 20 cm radio con- 
tinuum data. 



4.2. Star formation rate and local Schmidt law 

The radio continuum emission of the disk of M51 re- 
flects the spatial distribution of the current star forma- 
tion rate. In normal galaxies, such as M51, most of the ra- 
dio c ontinuum at 20 cm is non-thermal synchrotron emis- 
sion llCondonlll99l . radiated by cosmic rays interact- 
ing with the magnetic fields of the interstellar medium. 



The cosmic rays in turn are emitted by supernova rem- 
nants. The remaining radio continuum emission is free- 
free bremsstrahlung emission from thermal electrons in 
Hii-rcgions. The radio continuum thus traces the current 
star formation rate of massive stars. 

The observed optically thin radio continuum emission 
is known to be well correlated with FIR dust contin- 
uum emission in a large variety of sources, including nor- 
mal disk galaxies. The tight correlation was confirmed by 
iHelou et al.l l |l 98.^ using IRAS data and is described by 
the parameter 



q = log 



FIR 



[3.7510 12 Wm" 



log(- 



[Jy] 



(i) 



The parameter q was found to be 2.3 w ith a dispersio n 
of ~ 0.2 dex in a wide variety of sources llCondonill99l . 
We used this relation to deriv e the FIR flux at 15" reso- 
lution at each position in M51. lMurphv et all l(2005h have 
recently confirmed that this correlation holds for a map of 
M51 at 70 (im using MlPS/Spitzer at 19" resolution. In the 
disk, they find a median and scatter of qjo = 1.94 ± 0.19. 
Note that the correlation is expected to break down at 
yet smaller scales of a fe w hundred parsecs around mas- 
sivc star- forming regions ( Boulanger fc PeraultHl988h . 

The FIR inte nsity reflects the current star forma- 
tion rate (see e.g. lYun et al1l200lf) . as young stars form 
deeply embedded in their parental molecular clouds, be- 
fore dispersing their environment by forming H II regions 
and by supernovae explosi ons. Following the argument of 
iThronson fc Telescol ( 198fih . the FIR luminosity is propor- 
tional to the star formation rate: 



/. in; — SFR tpiR L IM 



(2) 



where L/M is the luminosity-to-mass ratio of the 
young stellar clusters. Assuming a typical disruption time- 
scale of tpiR = 2 10 6 yr and a Salpeter initial mass func- 
tion (IMF), the current star formation rate is SFR= 
6.5 10~ 10 (Lfir/Lq) Mq yr -1 . In summary, the star for- 
mation rate per unit area is related to the flux density at 
20 cm via 



E 



SFR 



[M Q pc" 2 Gyr' 



1.5310 s 



[Jybeam : ] 



(3) 



The gas depletion or consumption time is defined as 



y /y 



SFR- 



This is the time which would be 
needed to convert the total gas content into stars assuming 
that the SFR is constant with time and that there is no gas 
infall or recycling via stellar winds. Note that the inverse 
of the gas depletion time is the star formation efficiency. 

The global star formation rate of M51 is 2.56 M© yr^ 1 
and the global gas depletion time 0.76 Gyr. The latter 
value corresponds to a global star formation efficiency 
(SFE) of 13% per 0.1 Gyr. Both values, the depletion time 
and the SFE, agree to within a fact or of 2 with the val- 
ues derived by IScoville et al.l (|200l|) fro m Hubble Space 
Telescope (HST) Ha images and CO data. lMisiriotis et alJ 
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( 2006J) have recently studied the distribution of the ISM 
in the Milky Way. They find a similar SFR of 2.7 M Q yr" 1 
but factor ~ 5 larger gas depletion time of 3.57 Gyr in the 
Milky Way. Typical normal spiral galaxies need less time 
than found in the Milky Way, but more time than found in 
M51, to consume all the gas into stars. iKennicutti <)l998h 
finds a median depiction time of 2.1 Gyr for his sample of 
61 normal disk galaxies which shows variations between 
0.2 Gyr in starburst galaxies like NGC 5169 and 12 Gyr in 
early- type spirals such as M31. 

The SFR peaks above 100M Q pc _2 Gyr _1 in the center 
(FigJU, indicating a nuclear starburst, and drops radially 
to values of 3 M Q pc _2 Gyr _1 at 12 k pc distance. No t e that 
the C2 cluster complex studied bv iBastian et alJ l)2005[) 
shows a local SFR of 26OOM pc- 2 Gyr- 1 , while the other 
complexes studied by these authors show moderate local 
rates of ~ 60 — 70M Q pc~ 2 Gyr -1 . Regions of more than 
~ lOOM pc" 2 Gyr" 1 (i .e. 0.1 M^pc'V -1 ) are classi- 
fied as starbursts (e.g. iKennicutn 1 19981) . The center of 
M51 harbo urs a Seyfert 2 AGN surrounded by a ~ 100 pc 
disk/torus feohno et all Il996t) of warm and dense gas 
^Matsushita et al]ll 998i l2004h . 

Figure^] shows the variation of the radially averaged 
star formation rate in comparison with the Hi, H2 , and 
total gas surface densities in M51. Ignoring the center, 
the radial drop of the star formation rate closely resem- 
bles the drop of the total gas surface density. Neither docs 
the SFR show the rather flat distribution of the atomic 
gas nor the much steeper drop of the molecular gas. In 
contrast, we would have anticipated a good correlation of 
the star formation rate with the molecular gas since star 
formati on is known to oc cur only in molecular clouds. 
Indeed, iHever et alJ l)2004l) find a strong correlation be- 
tween the star formation rat e and the molecular g as sur- 
face density in M33. Similarly. IWong fc Blitd IMPOST) report 
a much better correlation of Ssfr with Sh2 than with 
Eht in their sam ple of 6 molecule-rich spirals. However, 
iKennicuttl ( 19981) who studied 88 galaxies and found that 
the disk-averaged SFR is much better correlated with the 
disk-averaged H 1 surface densities than with the H2 sur- 
face densities. In M51, Hi is often found downstream of 
the CO arms indicating that the Hi clouds are the rem- 
nants of GMCs photo-dissociated by young massive stars 
(Sec. 13.2(1 . The star formation rate may thus regulate the 
surface density of the atomic gas and hence explain the 
observed correlation. 

Figurc[|j] shows the variation of SFR in M51 with gas 
surface densities together with lines of constant gas de- 
pletion time. The gas consumption time in M51 varies 
between 0.1 Gyr in the center where the gas surface den- 
sity and the SFR surface density are high and 1 Gyr 
at large radii where the gas surface density and SFR 
density are low (Figs . HIBjl . In contrast, the Milky Way 
i Misiriotis et al.ll2006l) and the 6 C O-bright spiral galax- 
ies studied bv lWong fc Blit3 l|2002l) exhibit gas depletion 
times which are larger by a factor of about 10. They rise 
from 1 Gyr in the centers to ~ 10 Gyr and slightly more 
in the outskirts. 



Figure^] also shows that the star formation rate is pro- 
portional to a power of the total gas surface density im 
M51 



^-^SFR *—>p 



(4) 



i.e. follows a Schmidt-law (ISchmidtHl959l) . with nor- 
malization A = 1.1±0.5 and the slope n = 1.4±0.6 where 
Ssfr i s m um ts of M Q pc~ 2 Gyr -1 and ^ 
M pc- 2 . 



gas 



in units of 



In TableEl we compare the slope found in M51 with 
the results obtained for other galaxies and the results from 
simulations. The slope found in M51 agrees with the global 
Schmidt-law, seen in a study of disk-averaged and 
yiopp of 61 normal and 36 starburst galaxies 1 Kennicuttl 
Il998h . Note however that the slope of the 61 normal galax- 
ies is much less well-defined. Depen ding on the fitt ing 
method it varies between 1.3 and 2.5 (|Kennicuttll 1 998[) . 

In contrast, the observed slopes of local Schmidt laws, 
describing radial averages of £ gas and Ssfr in individual 
galaxies, do not in general agree wit h the global value, but 
vary strongly between 1.2 and 3.3. IWong fc Bht3 l|2002h 
studied radial averages of X) gas and X^sfr m a sample of 
spiral galaxies and derive local Schmidt laws with slopes 
between 1.2 and 2.1, assuming that extin ction depends 
on gas column density. iBoissier et alJ l(2003h study 16 spi- 
ral galaxies to study the local star formation laws and 
find a slope of ~ 2. In a similar study. IHever et all (|2004t) 
sho ws that M33 exhibits a significantly steeper slope of 
3.3. IMisiriotis et al.ll|2006[) find a slope of 2.2 for the Milky 
Way. 

4.3. Gravitational stability 

Exceeding a critical gas density may lead to the forma- 
tion of clouds and possibly stars. The critical density may 
in turn be determined by gravitation. Analogous to the 
well known Jeans stability criterion, the influence of dif- 
ferential rotation of a homogeneous, thin disk can be de- 
scribed by a dispersi on relation for axially symetric dis- 
turbences <Klevl2004) : a 2 = k 2 +o\^k 2 -2itG £ \k\. Here, 
a is the oscillation frequency, er gas the gas velocity disper- 
sion, k the wavenumber, and k the epicyclic frequency 
k2 = TT ("S + %)■ For a flat rota tion curve, h 2 = 2fl 2 . 
The system becom es unstable when the Toomre parame- 
ter l|Toomrelll96^ 



n r^- ^it(R) _ K(R)a sas (R) 
^ as[K) - £ gas (i?) - 7rGE gas (i?) ' 



(5) 



describing the competition of pressure and rotation on the 
one side and gravitation on the other side, at a given ra- 
dius, is less than 1, i.e. when the gas surface density ex- 
ceeds the critical surface density. 

In the following, we will assume that the disk of 
M51 is almost unstable as predicted by disk models (e.g. 



8 



Schuster et al.: A complete 12 CO 2-1 map of M51 with HERA 





A 


n 


Reference 


Local Schmidt-laws: 










M51 


1.12 ± 0.50 


1.35 ± 0.61 


this 


paper 


M33 


3.5 ± 66 


6.6 ± U.07 


Hever et al. 


f2004) 


Milky Way 




2.18 ±0.20 


Misirioti^^^i 


f2006) 


7 CO-bright spiral galaxies 




1.7 ±0.3 


Wong & Blitz 


f2002) 


16 spiral galaxies 




~ 2 


Boissier et al. 


f2003) 


simulations 


0.25 ±0.16 


1.31 ±0.15 


Li et al. 


f2006) 


Global Schmidt-laws: 










97 normal and starburst galaxies 


0.25 ±0.07 


1.4 ±0.15 


JKejmicrr^ 


f!998^ 


7 CO-bright spiral galaxies 




1.7 


Wons: & Blitz 


f2002) 


simulations 


0.11 ±0.04 


1.56 ±0.09 


Li et al. 


f2006) 



Table 2. Schmidt Law X)sfr = ^Sgas observed in M51 and other galaxies and samples of galaxies in comparison 
with the Schmidth-Law deri ved simulations . 5^sEfi anc ^ are m un its of M Q pc -2 Gyr" 1 and J2 sas m un its of M Q pc~ 2 . 
The values of the sample of IWong fe Blit d l)2002j) hold for an extinction correction that depends on the gas column 
density. 



iLin fc Pringleil987l) . i.e. <5 g as = 1, to calculate the critical 
gas velocity dispersion cr cr jt necessary to stabilize the gas 
against gravitational collapse. We will then compare <7 cr ;t 
with the observed velocity dispersion of the molecular gas, 
(TcOj as function of the galacto-centric radius. 



4.3.1. Rotation curve 

The rotation curve of M51, needed to calculate the 
epicyclic freque ncy and the critical gas densi ty (Eq.JSJ), 
was derived by iGarcia-Burillo et"aT1 (|l993bla|) from the 
CO 2-1 position-velocity diagram along the major axis 
of M51, using the velocities at the peak intensities, and 
correcting for the inclination. The rotation curve (Fig.Ql 
rises steeply within 10" of the center indicating the pres- 
ence of a compact nuclear mass component in addition to 
the central bulge and then stays constant at ~ 200 kms -1 
out to a radius of 10 kpc. We assume here that the rotation 
curve stays almost constant further out to at least 12 kpc, 
ignoring any effect of the companion galaxy. The defini- 
tion of the rotation curve used here appears to be more 
appropriate for the face- on galaxy M51 than the defini- 
tion used bv lSofuel l|l996j) who derived the rotation curve 
from the terminal velocities at which the peak intensity 
has dropped to 20% of its value, which leads to very high 
velocities of upto 260 kms" 1 at radii between 3 and 9 kpc. 



4.3.2. Critical velocity dispersion 

Figured also shows the radially averaged critical velocity 
dispersion cr crit derived from the gas surface density and 
the rotation curve, assuming Q gas = 1 (Eq.[5j). The critical 
dispersion, necessary to stabilize the gas against gravita- 
tional collapse, lies between 1.7 and 6.8 kms -1 . The disper- 
sion peaks near R = 5.5 kpc where the H2 surface density 
is high (Fig.^J. It then declines slowly to ~ 2kms _1 at 
R = 12 kpc. 



4.3.3. Observed velocity dispersion 

The velocity dispersion of the molecular gas is estimated 
here from the equivalent 12 CO 2-1 line widths Av cq = 
J Tdv/Tpk via ctco = Ai> oq /(2\/2 In 2). Figured shows 
the radially averaged velocity dispersion cr co . ^ drops 
from ~ 28kms~ 1 in the center to ~ 6kms _1 at radii of 7 
to 9 kpc. Further out, it rises again to values of ~ 8 kms" 1 . 
The rise of the observed dispersion at radii of 10-13 kpc, is 
due to the increased line widths in the companion galaxy. 

Similar vertical velocity dispersions are derived from 
Hi observations of spiral galaxies. Typically, these vary 
radially from ~ 12 — 15kms _1 in the central parts to ~ 
4 — 6 kms -1 in the outer parts (see review of iDib et alJ 
l2006tl . 

The CO velocity dispersion observed in M51 exceeds 
the critical dispersion at all radii. The ratio between the 
observed and the critical dispersions drop from a factor of 
~ 5 to almost 1 for radii between 1 and 5 kpc. This may 
indicate that the gas is stabilized against collapse at all 
radii. 

However, the observed gas dispersions are averages 
over the 450 pc beam and along the elliptical annuli, and 
are broadened relative to the intrinsic dispersions due to 
systematic motions which do not contribute necessarily 
to local support of the gas against gravity. Moreover, the 
influence of stars is neglected. 

5. Discussion 

5.1. Toomre stability 

Several authors have argued that galactic disks are self- 
regulated through gravitational instabilities to h ave their 
Toom re Q parameter of the order of 1 (see e.g. ICombesI 
As soon as gas dissipation leads to E gas > S cr i t at 
a given radius, the disk becomes gravitationally unstable, 
leading to density waves transferring angular momentum, 
causing dissipation and heating of the gas and thus an in- 
crease of the velocity dispersion, stabilizing the disk again. 
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H [kpc] 

Fig. 7. Left: The black drawn line shows the 12 CO 2- 
1 gas dispersion ceo derived from the observed equiva- 
lent widths. Positions without detection of 12 CO are ig- 
nored. The green drawn line shows the gas dispersion 
c C rit; assuming gravitational stability, i.e. S C rit = S gas . 
Ri ght: The dashed line gives the rotation curve derived 
bv lGarcia-Burillo et all (|l993a|) . 



Models indee d indicate that the d isk settles at the border 
of instability l)Lin fc Pringlelll987(l . 

However, in general the velocity dispersion and surface 
density of star s may play an i mportant role in the sta- 
bility analysis. iBottemal l|l993f) discuss the dispersion of 
stellar disks and the self-regulating mechanism that keeps 
Q near 1 over the entire disk of spiral galaxies. To take 
this into account, a total Q parameter has been used by 
e.g.lCombesI l(200l|): Q^t = Qgai + Qstar- Neglecting stars 
is valid only when the gas surface density divided by the 
gas dispersion dominates all other £/er terms. 

Observationally, Q gas has oft en been found to be 
larger than 1 by factors of a few. IKennicuttl l)l989j) de- 
rived radial averages of the Toomre-parameter of a sam- 
ple of 15 galaxies, assuming a constant velocity disper- 
sion of 6kms , and found that star formation is re- 
stricted to regions where Q gas < 1.5. A similar conclu- 
sion w as drawn from a study of 32 nearby spiral galax- 
ies by IMartin fc Kenni cuttl ll200l[) with few exceptions. 2 
Similarly. IWong fe"l3htJi)2002() also assume a constant ve- 
locity dispersion and find S gas w E cr it over a wider range 
of radii for some of their galaxies while others show devi- 
ations of factors of ~ 2. 

While in the inner 4 kpc of M51, the gas disper- 
sion is large and Qgas lies between 2 and 5, the stars 
may contribute here, and make the disk marginally sta- 
ble. Streaming motions are a consequence of spiral waves 



and gravitational instabilities l|Combesll200l[) . M51 shows 
strong streaming motions of 60 — 150 kms -1 in the plane 
of the galaxy in the inner spira l arms as can be s een in 
the interferometric CO maps of lAalto et al ] <|l999h . The 
broadening of the velocity dispersions due to streaming 
motions is naturally taken into account here. However, the 
observed dispersion is also broadened by the very steep, 
unresolved rotation curve in the center of M51. A more 
careful stability analysis would need to correct the ob- 
served dispersions for rotation gradients, take into account 
the influence of stars, and study individual regions of M51 
at high spatial resolution. 

5.2. The local Schmidt law and star formation 
efficiency 

Recent smoothed particle hydrodynamics (SPH) simula- 
tions of gravitational instability of isolated disk galaxies 
comprising a dark matter halo, a disk of stars, and isother- 
mal gas bv lLi et all l|2005Ll2006j) find a slope of 1.31 of the 
local Schmidt law (TablcEJ) , close to the slope of the global 
Schmidt law, which is also reproduced, and close to the 
slope of 1.35 found in M51. On the other hand, the slope 
of 1.3 is at the low-end of the local Schmidt law slopes 
observed so far. 

We find a normalization factor A of the local Schmidt 
law in M51 of 1.12 M w pc~ 2 Gyr~| (Table©, similar to A of 
gas-rich models in iLi et al.N2006h. This is a factor 5 larger 
than the average value of the lLi et alJ l)2006|) models. It is 
also a factor 5 larger than the no rmalization found in the 



large sample of lKennicuttl [ 

Interestingly, the normal galaxies of the IKennicuttl 



2 Both. lKennicutd il98Sft and lMartin fc Kennicuttl 11200111 in- 
cluded M51 in their sample, using FCRAO 50" CO 1-0 data 
combined with Ha data. However, the Toomre parameter and 
its variation with radius in M51 are not explicitly discussed. 



19981) sample show on average a factor ~ 3 longer gas 
consumption times than found in M51. In Sectio nLL2J we 
also noticed that the Milky Way JMisiriotis et alJl2006l) 
and t he 6 CO-bright galaxies studied bv IWongfcBlita 
()2002h exhibit consumption times which are larger than 
the consumption time found in M51 by about one order 
of magnitude. Apparently, the star formation efficiency of 
M51 is higher than in many other normal galaxies. 

The reason may be the interaction of M51 with its 
neighbouring galaxy NGC 5195, as has been proposed be- 
fore a s mentioned in the Int r oduction (e.g. iNikola et all 
120011 iHoward fc Bvrdl Il990l iToomre fc Toomrel 119721)" 
This i s supported by the recent simulations of ll*i et all 
()2004l) which show that interaction and merging can lead 
to instable disks and strong starbursts. 

6. Summary 

We mapped completely the interacting spiral galaxy M51 
in the 12 CO 2-1 line using HERA at the IRAM-30m tele- 
scope. The map includes the companion galaxy as well as 
the south-western arm out to radii of 12 kpc at linear res- 
olutions of 450pc (11"). These data were combined with 
maps of Hi and the radio continuum at 20cm at similar 
resolutions to study radially averaged surface densities. 
The star formation rate per unit area was estimated from 
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the radio continuum, allowing to study star formation laws 
like the Schmidt law. The critical density for gravitational 
instability was compared with the total gas surface den- 
sity along the radial averages. In detail, we obtain the 
following results: 

— Global properties: The total H 2 mass of M51 is 
1.941O 9 M and the global Hi/H 2 mass ratio is 1.36. 
The global star formation rate is 2.56 M Q yr _1 and the 
global gas depletion time is 0.8 Gyr. 

— The total gas surface density drops from ~7OM pc~ 2 
in the center to ~ 3M Q pc -2 at radii of 12kpc. The 
ratio of atomic to molecular gas surface density rises 
from 0.1 in the center to 20 in the outer regions. It is 

1 at ~ 4kpc. Up to ~ 6kpc the ratio follows a simple 
power law of R similar to that found in other spirals. 
At larger radii, the simple power law relation breaks 
down. 

— The gas depletion times in M51, vary between 0.1 
and 1 Gyrs, and are shorter than in many other non- 
interacting normal galaxies. Simulations indeed show 
that interaction can lead to high star formation effi- 
ciencies. 

— The star formation rate per unit area drops from 
~ 400 Mq pc~ 2 Gyr -1 in the starburst center to ~ 

2 M Q pc -2 Gyr -1 in the outskirts. It is much better 
correlated with the total gas surface density than with 
the surface densities of the molecular or atomic gas. 
The correlation follows a Schmidt law Esfr oc £™ as 
with an index of n = 1.4 ± 0.6. Only few studies of 
local Schmidt laws exist to date. The slopes of local 
Schmidt laws observed in other spirals vary strongly 
between 1.2 and 3.3. The slope of 1.4 agrees with re- 
cent S PH simulations of isolated disk galaxies lLi et alJ 
ll2005h . 

— The radially averaged dispersion of the observed 12 CO 
2-1 lines varies between 28kms _1 in the center and 
~6kms _1 in the outskirts. 

— We investigated the gravitational stability using the 
Toomre criterion. The critical gas velocity dispersions 
needed to stabilize the gas against gravitational col- 
lapse in the differentially rotating disk of M51, vary 
with radius between 1.7 and 6.8 kms . Observed ra- 
dially averaged dispersions derived from the CO data 
exceed the critical dispersions by factors Q gas of 1 to 5. 
Unresolved gradients of the rotation curve lead to an 
overestimate of the intrinsic velocity dispersion. The 
gravitational potential of stars not considered here, 
may lead to a marginally stable disk. 

In a forthcoming paper we will discuss the distribution 
of molecular gas, local Schmidt laws, and local critical 
gas surface densities, discerning the conditions in different 
parts of M51, i.e. the arm, interarm, and central regions. 
We will also study in more detail the spatial and kinematic 
distribution of atomic and molecular clouds, especially in 
the outer spiral arms. 
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